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Abstract 



, The lowest order QED radiative corrections to the radiative large angle Bhabha 

scattering process in the region where all kinematic invariants are large compared to 
the electron mass are considered. We show that the leading logarithmic corrections 
do not factorize before the Born cross section, contrary to the picture assumed in the 
J> I renormalization group approach. The leading and non leading contributions for typical 

' kinematics of the hard process at the energy of the $ factory are estimated. 
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O 1 Introduction 

00 

^ I The large angle Bhabha scattering process (LABS) plays an important role in e~^e collid- 
p || ing beam physics First, it is traditionally used for calibration, because it has a large 
p • cross section and can be recognized easily. Second, it might provide essential background 
1^ . information in a study of quarkonia physics. The result obtained below can also be used to 
construct Monte Carlo event generators for Bhabha scattering processes. 

In our previous papers we considered the following contributions to the large angle 
^ ' Bhabha cross section: pair production (virtual, soft 0, and hard 0) and two hard pho- 
c3 ■ tons 0. This paper is devoted to the calculation of radiative corrections (RC) to a single 
hard-photon emission process. We consider the kinematics essentially of type 2 ^ 3, in 
which all possible scalar products of 4 momenta of external particles are large compared to 
the electron mass squared. 

Considering virtual corrections, we identify gauge invariant sets of Feynman diagrams 
(FD). Loop corrections associated with emission and absorption of virtual photons by the 
same fermionic line are called as Glass-type (G) corrections. The case in which a loop 
involves exchange of two virtual photons between different fermionic lines is called Box-type 
(B) FD. The third class includes the vertex function and vacuum polarization contributions 
(rn-type). We see explicitly that all terms that contain the square of large logarithms 
ln(s/m^), as well as those that contain the infrared singularity parameter (fictitious photon 
mass A), cancel out in the total sum, where the emission of an additional soft photon is also 
considered. 

We note here that the part of the general result associated with scattering-type diagrams 
(see Fig. 1 (1,5)) was used to describe radiative deep inelastic scattering (DIS) with RC taken 
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into account in Ref. (we labeled it the Compton tensor with heavy photon). A Similar 
set of FD can be used to describe the annihilation channel . 

The problem of virtual RC calculations at the 1 loop level is cumbersome for the process 

e+(p2) + e'(pi) — . e+(p;) + e'{p^) + 7(A;i). (1) 

Specifically, if at the Born level we need to consider eight FD, then at the 1 loop level we have 
as many as 72. Furthermore, performing loop momentum integration, we introduce scalar, 
vector, and tensor integrals up to the third rank with 2,3,4, and 5 denominators (a set of 
relevant integrals is given in Appendices A, B). A high degree of topological symmetry of FD 
for a cross section can be exploited to calculate the matrix element squared. Using them, 
we can restrict ourselves to the consideration of interferences of the Born-level amplitudes 
(Fig. 1 (1-4)) with those that contain 1 loop integrals (Fig. 1 (5-16)). Our calculation is 
simplified since we omit the electron mass m in evaluating the corresponding traces due to 
the kinematic region under consideration: 

S ~ Si ~ —ti ^ —t —U ~ —Ui ~ Xi 2 ~ Xl 2 ^ ^^^5 (2) 
S = 2piP2 , t = -2p2P2 , U = -2P1P2 Si = > 

ti = -2pip[ , Ml = -2p2p[ , xi,2 = '^kiPi,2 , x'1,2 = 2A;ip'i,2 > 
s + Si + t + ti + u + ui = 0, s + t + u = x'l, 
si+t + Mi = -xi, t + xi = ti + Xi. 

We found that some kind of local factorization took place both for the G and B type 
FD: the leading logarithmic contribution to the matrix element squared, summed over spin 
states, arising from interference of one of the four FD at the Born level (Fig. 1 (1-4)) with 
some 1 loop-corrected FD (Fig. 1 (5-16)), turns out to be proportional to the interference 
of the corresponding amplitudes at the Born level. The latter has the form 

E, = {Ana)-'Y.\Mi\' = -f ^Mp'^OwPidw) ■ ^Tr(p27Ap) 

^\{u' + ul + s' + sl), (3) 



tXiXi 

" ' S Si U Ui 



X {u^ + ul + + si), 



tti\XiX2 X1X2 X1X2 X2X1 
4 r 4uiX2 



(47ra)-=^^Mi(M3* + M;) = -(1 + Z) — | 



Xi 



Au(si + ti)(s + t) 2 ,^ /X/ 

-|- ; \2sUUi + [U + Ui)(UUi + SSi — tti)\ 

X2Xi X1X2 

2 

+ -Y2tiuui + (n + Ui){uui + tti — ssi)] 

XiXi 

p[ + ki Pi-k ~ 

Oil' = 7p ) — 7m - 7/. 7p, Oiv = Ou'ip ^ /i), 

Xi Xi 
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where the Z-operator acts as follows: 



Z 



Pi < — *• s 

P2 < ^ P'2 U ^ 



It can be shown that the total matrix element squared, summed over spin states, can be 
obtained using symmetry properties realized by means of the permutation operations: 



{AixafF, F ={l + P + Q + R)^ = 
ssi(s2 + si) + ttiif + tf) + uuiiu^ + ul) 



(4) 
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X 



+ 

X1X2 

Eo + Oo- 



XiX'2 
-/n. 



SSitti 

t ti 



+ 



u 



+ 



Ui 



X2X2 XiXi X1X2 X2X1 



The explicit form of the P, Q, R operators is 



Pi < — ^ 

P2 < ^ 

ki ki 



-P2 
-P'l 



Q 



R 



P2 



-p'l 



P'2 P'2 



S < > Si 

t < — ^ tl 

U, Ui — > U, Ui 

s < — > tl 
Si < — > t 



pi,ki^pi,ki u,ui^u,ui 
Pi < — > —p'2 s < — > t 
p'l p'l si < — > tl 

P2,ki^P2,ki U,Ui^U,Ui 



(5) 



The differential cross section at the Born level in the case of large angle kinematics (|^) was 
found in Ref. M: 



do-o(pi,P2) 



a 



^ dVidVad'fci ' ' u^ 

F — 6^ > [pi +P2-Pi-P2-ki] 



^182^1 



(6) 



where ei, £2, and uji are the energies of the outgoing fermions and photon, respectively. 
The collinear kinematic regions (real photon emitted in the direction of one of the charged 
particles) corresponding to the case in which one of the invariants Xii x'i is of order yields 
the main contribution to the total cross section. These require separate investigation, and 
will be considered elsewhere. 

Our paper is organized as follows. In Sec. 2 we consider the contribution due to the 
set of FD Fig. 1 (5-8) called glasses here (G type diagrams). Using crossing symmetry, we 
construct the whole G type contribution from the gauge-invariant set of FD in Fig. 1 (5). 
Moreover, only the set of FD depicted in Fig. 2 (d) can be considered in practical calculations, 
due to an additional mirror symmetry in the diagrams of Fig. 2 (d, e). We therefore start by 
checking the gauge invariance of the Compton tensor described by the FD of Fig. 2 (d, e) for 
all fermions and one of the photons on the mass shell. In Sec. 3 we consider the contribution 
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of amplitudes containing vertex functions and the virtual photon polarization operator shown 
in Fig. 1 (13-15) and Fig. 2 (f, d). In Sec. 4 we take into account the contribution of FD with 
virtual two-photon exchange, shown in Fig. 1 (9-12), called boxes here (B type diagrams). 
Again, using the crossing symmetry of FD, we show how to use only the FD of Fig. 1 (9) 
in calculations. We show that the terms containing infrared singularities, as well as thoes 
containing large logarithms, can be written in simple form, related to certain contributions 
to the radiative Bhabha cross section in the Born approximation (|^). We also control terms 
in the matrix element squared that do not contain large logarithms and are infrared-finite. 
Thus our considerations permit us to calculate the cross section in the kinematic region (^), 
in principle, to power-law accuracy, i.e., neglecting terms that are 
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as compared to C(l) terms calculated in this paper. Note that the terms in (0) are less 
than 10-^ for typical moderately high energy colliders (DA$NE, VEPP-2M, BEPS). Unfor- 
tunately, the non leading terms are too complicated to be presented analytically, so we have 
estimation them numerically. In Sec. 5 we consider emission of an additional soft photon in 
our radiative Bhabha process. To conclude, we note that the expression for the total cor- 
rection, taking into account virtual and real soft photon emission in the leading logarithmic 
approximation, has a very elegant and handy form, although it differs from what one might 
expect in the approach based on renormalization group ideas. Besides analytic expressions, 
we also give numerical values, along with the non leading terms for a few points under typical 
experimental conditions. 

2 Contribution of G type diagrams 

We begin by explicitly checking the gauge invariance of the tensor 

u{p[)RZMPi)- (8) 

This was done indirectly in Ref. [^], where the Compton tensor for a heavy photon was 
written in terms of explicitly gauge invariant tensor structures. We use the expression 

d^fc jlxiP'i - k)'JaiPi -h- k)l\{Pi - hhf. 



_^ 7a (Pi - A;)7a(Pi -h- fc)7/.(Pi - ^)7a ] ^^^-^ 



ivr^ [ -Xi(0)(2)(g) 
1 - fci - fc) 
(0)(l)(2)(g) 



where 



{q) = {p,-k,-kf-m\ A2 = ^(L,,-i), i^xi=ln^. (11) 
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The quantity i?^^ corresponds to the FD depicted in Fig. 2 (d), while corresponds 
to the FD in Fig. 2 (e). The first term on the right-hand side of Eq. ([T0|) corresponds to the 
first two FD of Fig. 2 (d) under conditions @. The gauge invariance condition Rl1,k^ = is 
clearly satisfied. The gauge invariance condition regarding the heavy photon Lorentz index 
provides some check of the loop momentum integrals, which can be found in Appendix A: 

uip[)RZ,u{Pi)qae,{h) = Ak'^e^ih), A = -2^^^^^ - 6^^^^ . (12) 

Xi Xi 

The gauge invariance thus satisfied due to the Lorentz condition for the on shell photon, 
e(fci)fci = 0. As stated above, the use of crossing symmetries of amplitudes permits us to 
consider only R^^. For interference of amplitudes at the Born level (see Fig. 1 (1-4) and 
Fig. 1 (5-8)), we obtain in terms of the replacement operators 

iA\M\')a = 2W(1 + P + Q + m + Z)[E^^ + O^,^ - ^ - /f^l, (13) 

with 

E^i = }^hTr(p[RX^p,0u')-lTT{p,^,P2la), 

0^1 = ^ ^Tr(p;i?>^^pi7p) ■ 7Tr(M.p'2022'), (14) 

til 4 4 

m = ^ ^Tr(p;i?>^ipiOi2P27.P27p), 

= ^ ^Tr(p;i?>^ipi7pP27.P20r2'), 

^ p'l + ki pi~ ki 

Oiv = Ip ; — 7m -7m 7p, 

Xi Xi 

-P2 - ^1 -P2 + h 

O22' = 7m } 7p - 7p 7m> 

X2 X2 

^ Pi - ki -p2 + ki 
O12 = -7m 7p - 7p 7m' 

Xl X2 

p'l + ki -P2 - h 
Oi'2' = 7p ] — 7m + 7m ; 7p- 

Xl X2 

In the logarithmic approximation, the G type amplitude contribution to the cross section 
has the form 



dac = ^-{l + P + Q + m 



13 A 
--Lf + -Lu +2Lu In — 
2 *i 2 ' ' m 



(15) 
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3 Vacuum polarization and vertex insertion contribu- 
tions 

Let us examine a set of Fll-type FD. The contribution of the Dirac form factor of fermions 
and vacuum polarization (see Fig. 3) can be parametrized as (1 + Ff)/(1 — lit), while the 
contribution of the Pauli form factor is proportional to the fermion mass, and is omitted 
here. We obtain 

darn = ^ -2(1 + P + Q + i?)(Fi + n,)$, (16) 



F TT 



where 



r..^{(i„^-,)a-.,)-l.,-l.?.ic.}, (IT) 

In realistic calculations, the vacuum polarization due to hadrons and muons can be taken 
into account in a very simple fashion 0], just by adding it to Ut. 



4 Contribution of the B type set of Feynman diagrams 

A procedure resembling the one used in the previous section, applied to the B type set of 
FD (Fig. 1 (9-12a)), enables us to use only certain 1-loop diagrams in practical calculations, 
specifically three of those in the scattering channel with uncrossed exchanged photon legs: 



where 



and 



2^aV2 Re{l + P + Q + R) 



[(l-i^22')/g 



22' 



P2 * 
Pi 



-P'2 



P'vh 



s 

Si ^ 



u 



(19) 



JXI 



19 



TXi 
-'29 



X 



1 



16 1 



i7r2 (0)(g)((p2 + A;)2-m2) t 4 



X ^Tr {p2la{-P2 - k)jxP2-fp), 
d^k 1 



X 16 1 

-Tr {p2la{-P2 - khxp'2022'), 



d'k 



(0)(g) 



J 



4 1 



|--Tr(p'27,p;i?»piOi2^2(A + B 



(20) 



6 



7a(-P2 - k)jx 



B 



7a(-P2 + 

{—P2 + ky — m? 



Here 



+ 



+ 



7a(Pi -ki- k)-f^{pi - fc i)7^ , Jxipi -ki- k)jf,{pi - k)j^ 
-Xi{d) 
+ ki)l\{Pi - A:)7a 



{d) 



X'i(l) 
(pi - A:i - kf - 



m 



(q)^(p2-p'2 + kr-X', 
2 (1) = (p, _ ky - m\ (0) = fc" 



(21) 



Analytic evaluations divulge a lack of both double logarithmic (~ Li) and infrared logarith- 
mic (~ ln(A/m)L) terms in the box contribution. In spite of the explicit proportionality of 
the individual contributions to the structures Eq, Oq, and Jq, the overall expression terns 
out to be somewhat convoluted, despite its having a factorized form in each gauge-invariant 
subset of diagrams. We parametrize the correction coming from the B type FD as follows: 



1 1 <^ r A 

dcTfi = dcTo-L^AB, 

TT 



» „ , SS^ 2 , ^ ^ , , tt^ 

AB = 21n— + - + In— . 

uui h ssi 



(22) 



The total virtual correction to the cross section has the form 



da 



virt 



daa + d(7rn + dcr 



a 

TT L 



-Li 



+ L/^+41n- + AG + Arn + AB) +0(1) 
V 3 m / . 



(23) 



FV 



1/^2 ^2 ^2 ^2 

Ae + Arn = — l^ln h$i?ln h $q In — + $p In 77!- 

tti 



SSi 



f2 

h_ 

SSi 



tti 



where $p = P$, $q = Q^, and = R^. 



5 Contribution from additional soft photon emission 

Consider now radiative Bhabha scattering accompanied by emission of an additional soft 
photon in the center of mass reference frame. By soft we mean that its energy does not exceed 
some small quantity Ae, compared to the energy e of the initial beams. The corresponding 
cross section has the form 

dao ■ S'°^\ (24) 
_47ra r d^/c2 / pi ^ Pi ^ P2 p'2 V 

16tT^ J UJ2 V Plk2 Plk2 P2k2 ^2^2/ ^„<A£ 



d(7^°" = 
^soft 
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The soft photon energy does not exceed As <^ si = €2 = £ ^ £i ^ 62- In order to 
calculate the right-hand side of Eq. (El), we use the master equation M: 



Ana f d^k {qi 



Ana f (fk 2qiq2 



a , / Ae ■ m 
■ — m 

TT 



UJ 



IGtt^ J UJ {kqi){kq2 



j<Ae 



-iln^^"^ 
2 



£2 



TT 



h Li2 

3 



a 

TT 
COS" 



\ A ■ ei J 

( m\Aer 



Here we used the notation 



-hi 



(25) 



(26) 



In- 



gi,2 = (£1,2,91,2) 



-(91 - 92)^ > 



(27) 



where ei, £2, and 6' are the energies and angle between the 3 momenta qi,q2, respectively, 
and A is the fictitious photon mass (all defined in the center of mass system). 
The contributions of each possible term on the right-hand side of Eq. (El) are 



TT 



soft 



a 



Ai 

A12 

Ai/2' 

All' 
A22' 
A1.2 

Ai2' 



-Ai - A2 - a; - a;, + A12 + Ai/2' + All/ + A22' 

Al2' — Ai/2 , 

, Ae ■ m ^ , , Ae ■ m ^ , , Ae ■ m 
A2 = In ^ — , A[ = In — , A'o = In 



eX 

2L,ln ^ — + -L^ 

e\ 2 " 3 



e'lX 



4a 



Ls^ In 
Lti In 



/(A£- 

I / / 



V £^i4A2 
f{Ae-mf' 
\ e[eX^ , 
\Ae ■ raf\ 
ee'oX^ J 



2 2 [e. 



TT 



2 9i'2' 



, Li2 cos 
3 V 2 



TT 



+ Li2 COS' 



In' 



- y + Li2 I sin 



2 '^_2^ 

2 



^ , f{Ae-mf\ 1 2 1 2fe'A vr^ , . / 2 ^^i' 

i^i^ j + 2^- - 2 (7 j - y + i^^^ Y 

^ (Ae ■m)2 



L„ In 
L„ = In 



ee'nX'^ 



m 



In- 



2 

-Ml 



Lijiz) = - 



dx 



X 



ln(l — x). 



(28) 



where the center of mass energies of the scattered electron and positron, respec- 

tively; 611 , 62' are their scattering angles (measured from the initial electron momentum 
direction); and 6'i'2' is the angle between the scattered electron and positron momenta. 
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Separating out large logarithms, we obtain 

= ^^L. - 1) In + + + In + 

TT I AS UUi I J 

Cl'2' = COS ^i'2'- 



(29) 



This can be written in another form, using experimentally measurable quantities, the relative 
energies of the scattered leptons and the scattering angles: 



t 



1/1 + ^2 - 1 



— , Q=C0S6'-, -(l-Ci/2' 

£ 2 ■ yiy2 

1 + C2 U 1 - C2 tl 1 - Ci 



y2- 



2 ' s ^^'2 



-7=y'- 2 ' 



(30) 



— = z/1 + 1/2-1, — = yi^ — ■ 

s s 2 



6 Conclusions 

The double logarithmic terms of type and those proportional to Lsln(A/m) cancel in the 
overall sum with the corresponding terms from the soft photon contribution (^91). Omitting 
vacuum polarization, we obtain in the logarithmic approximation 



da 



soft+virt 



do-Q- 

TT 



L,(41n^ + Ai) + A(?/i,?/2,Ci,C2) 



(31) 



A, 



3 + lnil^l^ilil^ + ln.^^ + ^^ 



;i + ci)(i + c2) 



+ 



1/11/2 



t 

$ In h $p In — + <I>Q In — + In 

tti tti SS\ SS\ 



+ 2 In 



SSi 
UUi 



, ^ ^ \ 1 ^^1 
+ - $Q + $i? In . 



The function A(yi, 7/2, Ci, C2) is quite complicated. To compare it with Aj;^, we give their 
numerical values (omitting vacuum polarization) for a certain set of points from physical 
regions (|3^ ) and ?/i + ?/2 > 1, -D > (see Table 1). Considering the kinematics typical of 
large angle inelastic Bhabha scattering, we show the lowest-order contribution previously 
obtained and the radiative corrections calculated in this work. 

After performing loop integration and shifting logarithms [Li = + Lj^), one can see 
that the terms containing infrared singularities and double logarithmic terms ~ L^, are 
associated with a factor equal to the corresponding Born contribution. This is true of all 
types of contributions. 

The phase volume 

d^Kd^p^d^ ' I 

dr = — 5^ '{P1+P2-P1-P2 - h) 
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u 


yi 


y2 


Cl 


C2 


Al 


A 


1 


0.36 


0.89 


-0.70 


-0.10 


10.70 


-24.53 


2 


0.59 


0.66 


0.29 


-0.06 


4.86 


-11.41 


3 


0.67 


0.67 


0.50 


0.30 


5.82 


-35.58 


4 


0.68 


0.65 


0.60 


-0.50 


4.10 


-10.45 



Table 1: Numerical estimates of A^ and A versus ?/i, ?/2, Ci, C2 



can be transformed in various ways 0. We introduce the variables (see Eq. (|30|)) 

Ci = cos e'i , e'^ = p^p[ , 0<yi<l, -1 < ci,2 < 1, (32) 



which parametrize the kinematics of the outgoing particles (these do not include a common 
degree of freedom, a rotation about the beam axis). The phase volume then takes the form 

7rsdt/idy2dcidc2 , ^^r^^n^ \\ 

dr = — =B(z/i + 1/2 - l)Q{D{yi, 2/2, Cl, C2)), 

2V^(yi,2/2,ci,C2) 

D{yi, y2, Cl, C2) = -cl-cl- 2ci/2'CiC2, (33) 
p'^ = 2(1 -c..,) 

The allowed region of integration is a triangle in the yi, j/2 plane and the interior of the ellipse 
D > in the ci, C2 plane. 

We now discuss the relation of our result to the renormalization group approach. The 
dependence on Ae/e in ( ^T]) disappears when one takes into account hard two-photon emis- 
sion. The leading contribution arises from the kinematics when the second hard photon is 
emitted close to the direction of motion of one of the incoming or outgoing particles: 

^^hard ^ T^L^ ^ (daoizpi,p2,p[,P2) + dao{pi,zp2,p[,P2))dz 
Zn I i — z \ / 

+ dcTo Pi,P2, — ,P2 + dao Pi,P2,Pi, — dz2 

I ~ Zi \ ^1 / 1—^2 \ Z2 J . 

Z=l-X2, Zi = — ■ , X2 = — . (34) 

yi + X2 e 

The fractional energy of the additional photon varies within the limits Ae/e < X2 = uj2/e < 1. 
This formula agrees with the Drell-Yan form of radiative Bhabha scattering (with switched- 
off vacuum polarization) 

do-(pi,p2,pi,P2) = dxidx2V{xi)V{x2)dao\xipi,X2P2,—, — ] 

J \ Zx Z2) 

X V{zx)V{z2)dzxdz2, (35) 
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where the non singlet structure functions T) are |T0| 



(36) 



In our calculations we see explicitly a factorization of the terms containing double loga- 
rithmic contributions and infrared single logarithmic ones, which arise from G and m type 
FD. To be precise, the corresponding contributions to the cross section have the structure of 
the Born cross section (^. But the above claim fails to be true for terms containing single 
logarithms. Hence, the Drell-Yan form (^) is not valid in this case, and the factorization 
theorem breaks down, because the mass singularities (large logarithms) do not factorize be- 
fore the Born structure. That is because of plenty of different type amplitudes and kinematic 
variables, which describe our process. The reason for the violation of a naive usage of factor- 



ization in the Drell-Yan form has presumably the same origin with that found in Ref. |TT 
where the authors claimed that it is necessary to study independently the renormalization 
group behavior of leading logarithms before different amplitudes of the same process. Note 
that in the e/i — > 6/17 reaction, which can easily be extracted from our results, factorization 
does take place. We also see from (^) that factorization will take place if all the logarithmic 
terms become equal, i.e., ln(si/m^) = ln(s/m^) = ... . The source for the violation of 
the factorization theorem, we found, might have a relation to some of those found in other 



problems [p!2 



Numerical estimates (see Table 1) for the $ factory energy range {^/s ^ 1 GeV) shows 
that the contribution of the non leading terms coming from virtual and soft real photon 
emission might reach 35%. Additional hard photon emission will also contribute to 
and A. To get an explicit form of that correction, one has to take into account a definite 
experimental setup. 

Obviously, an analogous phenomenon of the factorization theorem violation takes place 
in QCD in processes like qq qqg and qq qq'y. A consistent investigation of the latter 
processes, taking into account the phenomenon found, can give a certain correction to pre- 
dictions for large angle jet production and direct hard photon emission at proton-antiproton 
colliders. 
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Appendix A 

Loop integrals for G— type Feynman diagrams 
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In this and the following appendices we used partially the results of our previous work |T3 



and refer to it for further details. After this digression let us turn to the problem. Two types 
of FD require different approaches. For the set of FD, labeled as glasses (G), only three 
independent external momenta are relevant due to the conservation law: pi + q = p'l + ki. 
Choosing pi, p[, q we use the notation: 

r (fk 1 _ r d^k 1 

J^... = J = gl^.y" + 4...P1P1 + %..p'vX + cl,„q^^^ (A.l) 

= / T^ (o)(i)(2)(g) = ^^i(^^i)'''' + ^^-^sAT^ + ^M'-' 

+ iTuiP^rpt + K222pTpiPi + K,,,q^q^q^ + K^^ipUr' 

+ K^,,{pWr^ + Ki2,{pip\qr\ 

where the inverse propagators are 

(0) = k'- \\ (1) = (pi - kf - m\ 

(2) = {p\ - kf - m\ (q) = ip[ -q-kf- m\ (A.2) 

A is a fictitious photon mass. The symmetrized tensor structures are defined as follows: 

[pqY'^ = pi^q^+fqi^^ (p^Y"^ = pf'p'q^ + pf'p^'' + p'p^", 
[gpY^P = gf'^pP + g^'Pp^ + gpyf", 
{pqrY"^ = p^q^r^ + p^q^" + p^q^r^ + p^q^^ + p^q^'r" + p^q^r". 

The vector and tensor integrals can be calculated by multiplying both sides of the ex- 
pression ( [A.2| ) by vectors Pi, Pi and q'^. Then one has to use the relations 

2p,k = (0) - (1), 2k,k = (q) - (1) + xi, 2p[k = (0) - (2), (A.3) 

and compare the coefficients before vector components on both sides. 

Considering the vector and tensor integrals with three denominators, we use ultra-violet 
divergent integrals with two denominators. Using the Feynman trick to join denominators, 
they can be expressed as 

d"fc 1 A2_ 

i7r2 [(A; -by-dy~^d 
d^fc kf" ,„/ A2 3 



i7r2 [(A; - 6)2 - \ d 2 



6^ In— - - . (A.4) 
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We put here the complete hst of these integrals (in approximation Eq.(0)): 

Joi = -^A + 1) Jlq = La — 1, J2q = La — Lt + 1, 

Joq = La — L^^ + 1, Ji2 = La — Lt^ + 1, J02 = -^a + 1, 



1 



01 



J 12 



L 



1 



Ji'q = K - 7;K) L 



2 
1 

2 

ip'i+pT)[^LA - ^L,, + ^ ), Jo"2 = pT{ ^La - 



lipt -k'i+ pT) {La -Lt + l), Jo'q = {P'l - {^La - \l^, + 



where 



Lq — Lt — In — Ly, 



1 ^1 T 



A2 



In — - ivr, La = In 



The scalar integrals with three denominators read 



- 2^ 



TT 



2L,L,, + L,^^ - - 



129 



2(x'i - Xi) 



02q 



Ji 



t + Xi 

1 



019 



-L 



L,{L, - L^J + -{U - L^,f + 2Li2 (1 + ^ 

-a;), La 



2X1 3x 



— , Li2(z) = -/— ln(l 



The coefficients for vector integrals with three denominators are 



O012 


— ^012 — 


O-Olq 


= Joiq + 


O02g 


= 0, bo 


0'12q 


t 

t-ti 


Cl2q 


tl 
t-ti 



—Lt^, C012 — 0, 
tl 
2 

— (L^, -1), 60] 
Xi 



1 ^ 

Xi 



Xi 2tLt (xi - t)L^, 

^02g + 7 , ,xo + 



Xi+t 



ixi + ty ixi + t) 



2 ' 



C02q 



Lxi ~ Lt 



Xi + t 



J I {t + ti)Lt,-2tLt 2 . _ ; ^ 



Jl2q + 



it + ti)Lt + 2tiLt, , 2 



The tensor integrals for G-type FD (see Eq.(A.l)) have the following form: 



T 

9oi2 



-012 



-(La - LfJ + - , aoi2 = O012 = 2f~^ " "^012 = ^ > 
/^012 = 7oi2 = O5 



T 



^Olq 



1 



(La-L^J + 



^Olq 



Lc 



Olq 



Xi 



Xi 



T 

-7oig 



2X: 



■(^xi-2), A 



T 
Olq 



-a 



9 
2^ 

T 

Olg " 



2X: 
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4(t + xi)"^^ 4(t + xi) ' ■ 8 
" 2(i + xi)' ^"'-^ " 2(i + xi)2 ^ 2(i + xi)' 



T 
9l2q 



I2q 

<T 
12q 

T 
ll2q 



0" 



1 J , tiLti — tLt 3 
T-^A — 7?:^ — n — r - , 



4{t -ti) 8 



Cl2q — 
T 



?,e + [ft I - ti 



3t^ 



At - ti 



y — "i; y — "i 
f + Atti + Stl^ , t{t - Ati^ ''^ ^2 



3ti t? 



Lu + 



-Lt + 



3t? 



- Atti - 3t? ^ Ati-t 
-Lt + — 



+ Atti + f\ t{t + 2ti) 2t + ti tti 

+ — TT^^t 



(A.9) 



2(t - 
ti(5t + ti) 



t{t + 5ti)^ , 3(t + ti) tti 



_ ti{t + 5ti) -t'^ + 5tti + 2tl t-lti 
- "TTTI T^-^ti ^71 T^n^ + 7^ 



2{t-hf 2(i-ti)3 
Four-denominator scalar integral reads: 

1 



2{t-hf {t-hf 



J 



12q- 



JoUq — 



-L,Lt, + 2Lt,L^, - Ll - 2Li2 f 1 - ^) - ^ 



6 



(A.IO) 



Vector 4-denominator integrals are: 
1 

CL012q — ^ 



&012, - ^ 
1 

C012, - ^ 



-(tXi + tiXl)Jl2q + (t + XlfJo2q " Xl(xi " ti)Joiq - ti{t + Xl)^ 
(^iXi + tXl)Jl2q - {ttl + x[Xl)J02q + Xl(Xl " tl)Joiq + tl(tl - Xl)Y 
-tlix'l + Xl)Jl2q + hit + Xl)^02q + Xl^Joiq " 



y = ^012 + XiJoi2q, d = -2tixix'i- 
2-rank 4-denominator tensors are: 



(A.ll) 
(A.12) 



T 
9oi2q 



]{Jl2q — XlCoi2g), 
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{t + Xlf{Jl2q - XlC012g) " (Xl^l + X'l^)«12g + Xl(^l " X'l)«0l9 



T _ 1 

- ^1 + Xi) (aoi2 + Xiaoi2g) 
1 



T 
l012q 



d 



{h - Xl) {Jl2q - XlC012q) + {Xlh + Xl^)&12g + Xl(Xl " ^l)Vg 



- (*1^ + XlXi)^02g + ti{ti - Xl)(«012 + Xl&012g) 
1 



d 



-tl{tl - Xl){Jl2q - 2X1C0129) + (X'l^l + Xl^)Cl29 - {ttl + XlX'l)c02g 



+ Xl(^l - Xl)&01g 
1 



012g 



d 



-{ttl + XiXi){Ji2q - X1C0129) + (x'l^i + Xi^)ai2g + Xi(xi " ^i)aoig 



+ ti{ti - Xl)(aoi2 + Xl«012g) 

1 



012g 



-0125 



d 



tliti + x'l){Jl2q - 2XlC012g) " (Xl^l + X'l^)Cl29 + (Xl + t) Co2q 



+ xi(xi -^OV? 

1 



t 



Jl2q - ^gm2q + ^l"oi2g + (x'l " ^l)/5oi2g + ^7oi2g 



We put now the coefficients of 3-rank tensor structures: 



(A.13) 



2g 



^[-(^ + Xi)'^i -ti{t + Xl)^8 + {th + XlXi)^18] 
\[{tti + XiXi)^i + - Xi)^ - (^1 - xO'^is] 



Kqg = 3Hl(^ + Xl)A-^M8 + ^l(^l-Xl)A8], 



'■95 

i^lll = 

i^ll2 = 

i^ll, = 

i^l25 = 

Klqq = 

K = 

^^qqq 

K122 — 



d 

^[-it + xifM - ti{t + xi)A9 + {ttl + xix'i)^i9], 

^[{tti + Xixi)^2 + ti{ti - xi)^9 - (ti - xO'^iq], 

+ Xl)^2 - tlAg + ti{ti - Xl)Al9], 



1 



t + Xi 

1 

t + x'i 

1 



[tiKiu + Oil^q - aoig - 2i^ig], 

[tlKiiq + - 



[^1-^lgg + Ci2q Cgig], 



* + Xi 

~l~[{tl — Xl)Kl2q + ai2g ~ "Olg 
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^299 = 




i^229 = 


tl 


-f^222 — 


tl 



T _ T ] 
"129 Co2gJ , 

T T 1 

7129 ^029] ) 

■'l2q %2q\ 1 



(A.14) 



where 



A T T A T T T A T T 

^1 = fl'129 ~ fl'029) ^18 = 5'012 ~ 5*019 + Xl9012q: ^8 = 5'l29 " fl'oi9> 

^2 = ^129 ~ ^i^lg; ^19 — ^012 ~ ^019 + X1Q0129) ^9 — ^129 ~ ^019- (A. 15) 

We give below some checking equations for coefficients before tensor structures of G-type 
integrals. The complete checking system can be obtained by contraction of general tensor 
expansion with relevant vectors, simplifying the numerators of the integrand and using a set 
of vector integrals given above. Additional check can be inferred by contraction with metric 
tensor. In this case the scalar integrals should be used. The complete set of 10 equations 
for the 2-rank tensor and 24 equations for the 3-rank 4-denominator tensor integrals for the 
G-type was convinced to be fulfilled. For definiteness we give four equations of such a type, 
obtained by contraction with metric tensor. They are: 

45J129 + ^^0129 - ^l"0129 + (Xl - tl)Pm2q + {t + Xl)lQl2q = Jl2q, 

6Kig-tiKn2 + {Xi~ii)Knq + tKiqg + {t + Xi)Ki2q = ai2q, 

6K2g -tlKi22 + iXl +t)K22q + tK2qq + iXl -tl)Ki2q = 6129, 

QKgg + tKgqq+{Xl-tl)Kiqq + {t + Xl)K2qq-tiKi2q = Cuq- (A.16) 

Another indirect check is the absence of infrared divergence containing terms in all the vector 
and tensor integrals. 

Appendix B 

Loop integrals for B— type Feynman diagrams 

We use here the following set of denominators: 

(1) = {pi-ky-m\ {2) = {p-,-k,-kf -m\ {Z)^{p2 + kf-m\ 

(4) = {pi-ki-p'^-kf-\\ {b)^k^-\\ (B.l) 

4-momentum conservation law we use reads pi +p2 = p[ +P2 + ki. Scalar products of the loop 
momentum k with the external 4- vectors can be expressed in terms of the denominators: 

2piA; = (5) - (1), 2p2A; = (3) - (5), 2^;^ = (4) - (2) - i - xi, 

2kik = (2) - (1) + xi, 2p'^k = (3) - (4) + t. (B.2) 
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Using these relations one can consider only one type of integrals with 5 denominators, namely 
the scalar one. Using the elegant technique developed in the paper of Van-Neerven and 
Vermasseren it can be expressed in the form: 



</l2345 

D2 



— [-Di J2345 



-D2t/l345 + -D3J1245 + -D4J1235 + -D5J1234] 
' , , , 1 



D = 2ssitxiXi, 



Sit[-t{s - Si) - sxi - SiXi - XiXij, 
st[t{s - si) + sxi + six'i - XiX'i], 
XiX'i[-tis + si) - sxi + Six'i + XiXi] 
sXiHs - Si) + sxi - Six'i - XiX'i], 
SiX'Ms - si) - sxi + Six'i + Xix'i]- 



(B.3) 



It is interesting to note that the method described above to calculate the coefficients of the 
tensor structures cannot be applied to the tensor integrals with 5 denominators given above. 
Some additional information is needed to close the system of algebraic equations. 

We mention a trick which permits to obtain additional equations for vector and tensor 
integrals whose denominators do not contain the term k"^ — It consists in shifting a loop 
momentum. Thus, for Jf234 we have 



k 



ivr^ (1)(2)(3)(4) 



k=p\—k 



d^k {pi - k) 

i^(i)(2)(3)(4) 



Pi -^1234 + a(Pi + P2) + cki + dp[ 



(i) = P-m2, {2) = {k-kif 



m 



(3) = ipi+P2 + kf - , (4) = {k-p\ - k^f 



The comparison of right hand side of this equation with the standard expansion 

0^1234 = {aPi + bp2 + cki + dp\)t2M 

leads to the new relation: 

^1234 = <^1234 = &1234 ■ 

Analogous useful relations can be obtained for tensor integrals as well. We put below the rel- 
evant scalar, vector and tensor integrals with 3 and 4 denominators from (B.l) and introduce 
the parametrization: 

d^k 1 rd^k k" 

[aij...pi + bij,„p2 



J i' 



ivr^ (i)(j) . . . 



'Jij... 



+ Cij..ki + dijp^)^, 
d^k k^'k^ 



{g'^g + a^'pipi + h^'p2P2 + c^'kiki + d^'p[p\ 



i7r2 . . . 

+ oP^ivm) + I3^{piki) + 7^(PiPi) + P^{PiP2) + (T'^{kiP2) + r^(p'i^i)) 
Vector 3-denominator integrals are: 

Lt 



:b.4) 



ij... 



^245 — — C245 - 

Xi 



J 2 



'245 



+ 



d 



245 



t + Xi 



-J' 



245 



+ 



&245 = 0, 

(Xi - h)Lt 



t + Xi 

(t+xir ' (t+xir 
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^145 = 0, Ci45 = du5 - 



Xi 



_ _ _ 2Lt 

Q345 — ~C345 — —0345 — — , O345 — — J345 H —, 

O125 — J125 -\ bi25 — di25 — 0, C125 = , 

Xl Xl 

O235 = ~C235 = , 0235 = 0, 

S-X2 

O235 — <>'235 — 7 + 7 \o-^si, 

S-X2 [s- X2Y {s - X2r 
Ls 

^135 — —bis5 — , C135 = (ii35 = 0, 

S 

•2234 — ~C234 — J234 ; O234 — ; "234 — ~<^234 H 

Si Si Si 



"^51 — 

O123 — J123 + bi23, bi23 — , "123 — 0, 

S — Si 

_ s ^ 2 2,sl,, (s + .Si)I. 

C123 — <^123 r 



S — Si S — Si (s — Si)2 (s — Si)2 ' 

■^x'l ~ 2 -^x'l 

O124 — Jl2A, bi24 — 0, C124 = —Jl24 H —, , 0^124 = r , 

Xl Xl 
s , , 2x[L^>-{s + x'i)Ls , _ J 

fll34 — 7<^134 -\ 7 7T5 , ^134 — "134 — J134, 

s-Xi {s-Xi) 

C134 — ai34 — 7 '^134 i 7 7T5 ■ V^-'^J 

s- Xi [s- XiY 



Vector integrals with 4 denominators read: 



01245 


Asa 
A3 


, OI245 — U, CI245 — A , 

^3 




C^1245 


Asd 

A3 


, A3 = 2iixix'i, 




A3a 


= X'l[Xl(2^1 + X'l)jl245 + x!lJl2A - XlJl2h -{t + Xl) J245 


+ {ti+Xl)Jl45, 


A3C 




XlX'l''^1245 + X'l>^124 + Xl''^125 " + Xl)>^245 + " 


■X'l)Ju5\, 


A3d 


= xi[- 


-XlX'l'^1245 - xUl24 + Xl'^125 + (Xl " ^l)^245 + {t 


-Xi)'^145]- 



(B.6) 
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A4a A46 


A4c 


Ol235 


— A ) 0i23b — A ) C1235 
ZA4 ZA4 


A4 ' 




= 0, A4 = 2sxiX2, 




A4a 


= X2[SX1^1235 - (S - Sl)^123 - 


(S - X2)^235 + Xl'^125 + 5^135], 


A46 


= Xl[sXl-^1235 + (S - Sl)-^123 - 


(S + X2) ^235 - Xl>^125 + S-^135], 


A4c 


= s[-sxi^i235 + (X2 - xO-^m 


+ (S - X2)-^235 + X1J125 - sJi35, 



A2a , A26 A2C . „ 

«1345 — — T— , O1345 — C1345 — 01345 — — — , ZI2 — ZStU, 

A2 A2 A2 

A2a = -St{s + t)Ji345 +t{s + t)J345 + S{S + t)Ji35 + {ut - Sx'i)Ju5 

+ {US - tXi)Jl3i, 

A26 = -St(s + m)Ji345 + ^(S - m)J345 + S(S + m)Ji35 - (S + m)Vi45 

+ («Xi - S^)-^134, 

A2C = s[stJi345 -tJs^Q - sJi35 + {S + u)Ji45 + {t - u)Ji3i]. (B.8) 



^2345 — — C2345 — — — , O2345 — ; "2345 — — , 

Z\i Z\i Z\i 

Ai = — 2si-uit, 

Ala = -SlMlt-^2345 - + Xl)'>^245 - Ml'SlJ234 + - X2)-^235 + ^MiJ345, 

Alb = -Slt(t + Xl)'>^2345 + (^ + Xl)^<^245 + Sl(^ + Xl)<^234 + (miXi + -Si^)<^235 
+t{ui - Si) J345, 

Ale = - X2)-^2345 + («lXl + Sl^)<^245 + Si(Mi - t)J234 + (S - X2)^-^235 

+^(s-X2)J345. (B.9) 



CH234 — -^1234 + -T— , &1234 — "T" , Ci234 — — </l234 T 1" -7—, (B.IO) 

(^1234 = -./1234 + -T^ - -r^, A5 = 2siXiX2) X2 = s - «i - X'l, 

^5 ^5 

Asa = X2[-{S - Sl)Jl23 + {s - x'l) Jl34 + Xi-^124 " S1J234 + SlXi-^1234], 
Asfe = X'llis - Si) J123 + (2Si - S + Xi) J134 - xi-^124 - S1J234 + SlX'iJl234], 
Asc = ■Si[(x2-Xi)</l23- (■S-X'l)</134 + Xi 0^124 + S1J234- SiX'i^l234]- 

We put now the tensor coefficients for B-type integrals with 4 denominators. 

fl'l245 = i [2»^124 - 0124 " XlCl245 + {t + Xl)c^l245] , 
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^1245 — T [x'l ("•^124 + 0124 — C145) + tiai45 — (t + Xi)a245 

+ ^lXiai245 - X'l{t + Xl)c^l245], 

1 

X 

1 



-1245 - T [hi-Jl24 + 0124) + X1C125 + (^1 - Xl)ci45 " XlX'lCi245] , 

XiXi 



^1245 = -1 7 [Xl(-'^124 + 0124 " ^245) + (*1 " Xl)ci45 " ^1^^245 " XlXWl245] 

flXl 

/^1245 — [—Jl24 + ai24 + C145 + XlCms] , 

Xl 

7l245 = 7- [■^124 - Ol24 + ^245 + Cl45 + {t + Xl)c^l245] , 

'^1245 — ~T ["-^124 + ^245 + Xl^mS ~ + Xl)c^l245] , 
Xl 

^1245 ~ *^1245 ~ P1245 ~ '^1245 ~ ^ • 

As a check one can use the result of contraction by the metric tensor: 

45'l245 + Xl/5m5 ~ ^l7l245 + X'l'''m5 = ■^124- 



fi'1235 — 2[^Jl2'i — «123 + ^123 — X1C1235] , 

«1235 = — [X2^123 - (Xl + X2)0l23 + Xiai25) " X1X2C1235], 
SXl 

^1235 = — [Xl(^123 - 0235) + (Xl + X2)&123 " X2&235 " X1C1235] , 

SX2 

C?235 = [■S('^123 + &I23) - (S - X2)a235 + X2C123 " SX1C1235] , 

X1X2 



1 



a 



1235 



— -[~>^123 + Ctl23 ~ Ct235 ~ ^123]) 



S 

1 



/^1235 ~ [~<^123 + CH23 + X1C1235]) 

Xl 

"^1235 = [— <^123 + ^235 — ^123 + X1C1235], 

X2 

^1235 ~ 7l235 ~ P1235 ~ '''1235 ~ ■ 

One of the checking relations here has the form 

4^^235 + S«?235 + Xl/5?235 + X2(7f235 = >^123- 
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5'l345 — 2^Jl34. + tCi345], 



n345 



[{S + t) Vi34 + t{Xl - S- t)ai45 - (s(s + t) + txi)a^l34 



stix'i - s-t) 

+ Xl{s + t){ci45 - C134) + t{s + tfci345], 
^1345 = -[&134 - &345 - (X'l - ^)Pm5]> 



^1345 



= O^Ls = n^45 = ^(^^7-3^— ^[(Xi - ^)(C145 - C134) - S(6l34 - tc^s)], 



"1345 = ^^^^rrj3^H(xi-S-t)a345 + x'l(x'l-^)(ci45-Ci34) 
- sxi(ai34 - -^134) + S^Xici345], 

= 7L5 = ^(^/^ -s-t) ^ ^^^^"^ ~ ^^^^^^-^ ~ '^'^ ^^^''^ ~ ^^^''^^ ' 
pLs = ^^Ls = ^^(^^rr^— ^[-(xi - ^)^ci45 + t{x[ -s- t)as45 

+ (x'i(xi -t)- st)cis4 + s{x[ - t)bis4 - st{x[ - i)ci345]- (B.15) 
The relation of the same type for the above coefficients reads: 

4^L5 + xicLs + ^"Ls + (xi - ^O/^Ls + (x2 - -"Of^Ls = -^m- (b.i6) 



5*2345 ~ 2 ["^234 + Xl'^2345 + + Xl)'^1345]) 

^2345 — ^2345 — ~/^2345 — ^["^^1345 — (si + Xl)a235 + •5l^Ct2345], 

b2345 = —n T -7t[si^(^235 " ^345) " Xl(Xl + t)a235 - t{t + Xl)«345 



^^2345 



- Sit(xi +t)&2345], 
1 

Xi + si + t 

+ t{xi + Si)a345 + (Xi + Si)^a235 - si^(xi + 51)02345) 



C?245 — C?234 TT^ — ; —7 ( Sit(a245 " ^234) 

sit{xi + si + t)^ 



•^2345 — ~^2345 — J^[~Xl0235 — ^«345] ; 

72345 = -'^2^45 = —Tf T [5l^(«245 " ^234) + ^(Xl + Sl)a345 

Si'CyXi + Si + rj 
+ (Xl + Sl)^«235 - sit{xi + ^1)02345] 

Plub = , -7y[-Si^«234 +Xi(Xi + Si)«235 + ^(Xi + Si)«345 

siHXi + si + t) 

- 51^X102345 - Sl^(Xl + ^)c^2345]- (B-17) 
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The above coefficients fiave to satisfy ttie reiation 

- Xl4345 + - X2)Qi2345 " + XOtLs " ^iPLs = ■^234- 



T 
fl'1234 



■^123 - Xl 



"■1234 
iT 

^1234 



--1234 



"1234 



A(2) 

2— T h J1234 + ^1234 , 



1234 , 

- 2- 
2- 



A A 

A(2) A(i) 



+ •^1234 + &1234 + C1234 — 271234 , 



•^1234 — 



/3; 



1234 



T 
7l234 



T 
P1234 



A 

A(2) 

— I- 01234 , 

A(=^) A(2) 

- 2 

- 2 



+ <^1234 + ^1234 + S1234 — 2q; 



1234 , 



A A 

A(i) A(2) 



A 

A(2) 

A(2) 



~ ■-^1234 ~ ^'1234 + 7l234 
— </l234 — &1234 + '5l234 



^^1234 — 



A 

^'1234 + <5;i234 , 
&1234 + 71234 , 



'1234 



A(2) A(i) A(3) ~ ~ . . 

2 — ^ ^ -. h ^1234 + &1234 + /3l234 ~ 0^1234 ~ 7l234; 



AAA 
wliere tfie quantities witfi tfie sign ~ are defined as follows: 
1 



(B.18) 



O1234 — — T\Ls 



^1234 — —r 

X2 



,A(2) 



Xi- 



+ 



X'l 



Xi 



J Jim H 



xi ^ si(s - xi) 



Xi(si -X2) 



C1234 — — r 
X2 



A(2) si , , si , si . 
X2 ^ X2 X2 s-si 



</l23 + 



2-L, 
s — Sl 



2-L 



xi 



2si 



xi 



(s-si)' 



-(Ls — LsJ 



ttl234 



si(s - xi) ' 



/3: 



A(3) 



1234 



71234 = — T 
X2 



Xl—. -'123 H 



+ 



Si 



Xi(s - si) 
s-Xi 



(B.19) 
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One of the checking relations takes the form 

25'l234 + Xl/^1234 + X2(^l2U + ^1^12Zi — ^134 — ^234 + Xiai234- 

At the end of this Appendix we give the table of scalar integrals with two, three and four 
denominators. We imply the real part everywhere and the ultraviolet asymptotic is assumed 
as well. 



where 



Jl2 


= -1 + , 


Ji3 = 1 + La — i^s , 


Jl4 


^1 + La 


- ^xi 


, Ji5 = J24 = ^34 = -'35 = -^A + 1 , 


■-'23 


= 1 + La 


~ Lsi 


, J25 — 1 -\- La — , 


Jib 


^1 + La 


-Lt, 




La 


= In—, 


Ls = 


In^, LA = ln— , 


Lsi 


= ln— , 




= ln— , L^i^ln— , Lt = ln— . 



3-denominator scalar integrals are 

1 



J 123 
Jl2i 

J235 
J135 
J245 
Jl45 



2(5 



-4-- 

2 ^1 6 



<^345 
•^125 = 



^L't + 



1 

Xi 



27r2 



TT 



2 ^1 3 





1 r 


s - 


-xi [ 




1 


Si +Xi 


1 

s 






1 n 


i + Xi 




1 



1 



^Li + -L^, - 2L,L^; + 2Li2 | 1 - ^ 

^ r 2 r 2 

2^s, + 2^X1 



27r2 



-2L,,L^^+2U2 (l + f^) 

, -^234 



-t + x'i 



1 

2^X1 
1 



-Ll-^Ll^ +2Li2(l + 



1 

Xi 



2 
27r2 



TT 



X'l 



-Lt/ - -L^ 2Li2 1 - — 

2 ^1 2 * 2 I t 



4-denominator scalar integrals read 

1 



J, 



1245 



XlXl 



1 

-'2345 = — 7 



-f'si ~ Ls^L\ — 2Ls^L^^ + 2LsjLt — ^ 



2L^^L( + 2Ly,'^Lt 



27r2 



(B.20) 
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(B.22) 



23 




1 2 

Ls~ LgLx - 2LsL^'^ + 2LsLt + 7— , 
^ Ll + LsL, - 2L,L^, + 2Li2 (l - ^) - 5^ , 
-Ll - + 2L.,Lx; - 2Li2 (l - -) - 7^ 



(B.23) 
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